function index=GetBoundarys(data)
[x,y]=find(data~=0);
k=boundary(x,y);
index=[x(k) y(k)];
% temp=zeros(60,60);
% for i=1:length(k)
%      temp(index(i,1),index(i,2))=data(index(i,1),index(i,2)); 
% end
% boundarys=temp;
end